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Wind braking of magnetars 

H. R. X. Xu2, L. M. Song^ and G. J. Qiao^ 

ABSTRACT 

Considering recent observations challenging the traditional magnetar model, 
we explore the wind braking of magnetars. There is evidence for strong multi- 
pole magnetic fields in active magnetars, but the dipole field inferred from spin 
down measurements may be strongly biased by a particle wind. Recent challeng- 
ing observations of magnetars may be explained naturally in the wind braking 
scenario: (1) The supernova energies of magnetars are of normal value; (2) The 
non-detection in Fermi observations of magnetars; (3) The problem posed by the 
low-magnetic field soft gamma-ray repeaters; (4) The relation between magne- 
tars and high magnetic field pulsars; (5) A decreasing period derivative during 
magnetar outbursts. Transient magnetars with < — -Erot may still be mag- 
netic dipole braking. This may explain why low luminosity magnetars are more 
likely to have radio emissions. A strong reduction of dipole magnetic field is 
possible only when the particle wind is very coUimated at the star surface. A 
small reduction of dipole magnetic field may result from detailed considerations 
of magnetar wind luminosity. In the wind braking scenario, magnetars are neu- 
tron stars with strong multipole field. For some sources, a strong dipole field may 
be no longer needed. A magnetism-powered pulsar wind nebula will be one of 
the consequences of wind braking. For a magnetism-powered pulsar wind nebula, 
we should see a correlation between the nebula luminosity and the magnetar lu- 
minosity. Under the wind braking scenario, a braking index smaller than three is 
expected. Future braking index measurement of a magnetar may tell us whether 
magnetars are wind braking or magnetic dipole braking. 
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1. Introduction 

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are magnetar 
candidates, i.e., neutron stars powered by strong magnetic field decay (Thompson & Duncan 
f995, f996). In studying them, the assumption of magnetic dipole braking is often employed 
(Duncan & Thompson 1992; Kouvehotou et al. 1998). However, the magnetic dipole braking 
mechanism is originally designed for rotation-powered pulsars. Since both the persistent and 
burst emissions of magnetars are from a different energy reservoir (magnetic energy instead 
of rotational energy), it is possible that they have a different braking mechanism, e.g., wind 
braking (Harding et al. 1999; Thompson et al. 2000). 

A strong dipole magnetic field obtained by assuming magnetic dipole braking is often 
taken as confirmation of a neutron star's magnetar nature (-Bdip > -Bqed = 4.4 x 10^^ G, 
Kouveloitou et al. 1998). However, the magnetic dipole braking assumption will also result 
in several problems challenging the magnetar model (Mereghetti 2008; Tong & Xu 2011). 

1. The spin down time scale of a newly born magnetar will be less than the shock breakout 
time due to the presence of a strong dipole magnetic field. This will cause the super- 
novae associated with magnetars more energetic than canonical supernovae (Duncan & 
Thompson 1992). However, observations of supernova remnants associated with AXPs 
and SGRs show that the corresponding supernova energies are of canonical value (Vink 
& Kuiper 2006). This failed prediction of the magnetar model may be circumvented 
if the initial rotational energy of magnetars are carried away in non-electromagnetic 
form, e.g., gravitational waves (Dall'Osso et al. 2009). However, in Dall'Osso et al. 
(2009), a relatively low dipole magnetic field is also required (-Bdip < 10^^ G). If mag- 
netars have a different braking mecahnism and consequently their dipole magnetic field 
is much lower, this may explain their supernova energy problem. 

2. If AXPs and SGRs are neutron stars with strong dipole field, then although they rotate 
rather slowly (periods: 2-12 seconds), they will also accelerate particles to very high 
energy. In the outer magnetosphere, these particles will emit high-energy gamma-rays 
which are detectable by Fermi-LAT (Cheng & Zhang 2001). This may be viewed as 
an independent measurement of strong dipole magnetic field, i.e., through unipolar 
induction effect. However, Fermi-LAT observations of all AXPs and SGRs show no 
significant detection (Sasmaz Mus & Gogus 2010; Abdo et al. 2010). Therefore, there 
are conflicts between the out gap model in the case of magnetars and Fermi-LAT 
observations (Tong et al. 2010a, 2011). It is possible that magnetars have a different 
braking mechanism and their dipole magnetic field is not so strong. 

3. In the traditional picture of the magnetar model, magnetars are young neutron stars 
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with both strong dipole field and strong multipole field (Thompson & Duncan 1995, 
1996; Thompson et al. 2002). The observation of the low magnetic field soft gamma-ray 
repeater SGR 0418+5729 has challenged the traditional magnetar prescription (Rea 
et al. 2010). This source tells us that magnetar-like activities (an anomalous X-ray 
luminosity or SGR-typc bursts) do not require a strong dipole magnetic field. The 
timing of SGR Swift J1822. 3-1606 further strengthens this point (Rea et al. 2012a). 
The originally required strong dipole magnetic field in most AXPs and SGRs mainly 
provides the braking torque. It is possible that not only SGR 0418-1-5729 but also 
many other AXPs and SGRs do not have a strong dipole magnetic field if they have a 
different braking mechanism. 

4. There are high magnetic field rotation-powered pulsars (HBPSRs) along with magne- 
tars (Ng & Kaspi 2011). Although they are close to each other on the P-P diagram, 
they show very different timing behavior. The timing behaviors of HBPSRs are similar 
to that of normal pulsars (Ng & Kaspi 2011). Therefore, it is reasonable that they 
have the same braking mechanism as that of normal pulsars. However, magnetars are 
very noisy (Gavriil & Kaspi 2002; Woods et al. 2002; Archibald et al. 2008), and the 
period derivatives of magnetars can vary significantly (up to a factor of 10, Gavriil & 
Kaspi 2004; Camilo et al. 2007; Woods et al. 2007). Therefore, it is possible that 
magnetars have a different braking mechanism, e.g. wind braking. The variation of 
wind luminosity will cause the variation of their period derivatives. 

All these issues are related to the dipole magnetic field and the braking mechanism of 
magnetars. A different braking mechanism of magnetars may help to solve these problems. 
Electrodynamics of magnetars show that they may have globally twisted magnetospheres 
(Thompson et al. 2002; Bcloborodov & Thompson 2007). The twisted magnetosphere 
will enhance their spin-down torque and also modify their persistent emissions. Changes 
in their global magnetosphcric structure will result in changes in their spin-down rate and 
persistent fiux (Bcloborodov 2009). In the case of wind braking, the large scale dipole field 
is unchanged. It is the change of the particle wind luminosity that causes change of the 
spin-down rate. Since both their persistent emissions and the particle wind are magnetism- 
powered, it is natural that their spin-down behavior and persistent emissions are correlated. 
Based on previous researches (Harding et al. 1999; Thompson et al. 2000), we explore the 
wind braking of magnetars in more detail and apply it to all AXPs and SGRs. A comparison 
with up-to-date observations is also given. 

Observations supporting the existence of a particle wind in magnetars are given in 
Section 2. Rotational energy loss rate due to a particle wind is calculated in Section 3. Several 
aspects of wind braking of magnetars are given in Section 4. Discussions and conclusions 
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are presented in Section 5 and 6, respectively. 

2. Existence of a particle wind 

2.1. Qualitative description of wind braking of magnetars 

In the magnetic dipole braking scenario of normal pulsars, the star's rotational energy 
is carried away by magnetic dipole radiation plus a rotation-powered particle wind (Michel 
1969; Xu & Qiao 2001; Spitkovsky 2006). The rotational energy loss rate is quantitatively 
similar to the magnetic dipole radiation in vacuum (Xu & Qiao 2001; Spitkovsky 2006). The 
surface dipole magnetic field is almost the same as that of magnetic dipole braking in vacuum. 
A particle wind mainly causes higher order modifications of pulsar timing, e.g. braking index 
(Michel 1969; Manchester et al. 1985; Xu & Qiao 2001; Contopoulos & Spitkovsky 2006; 
Wang et al. 2012) and timing noise (Lyne et al. 2010; Liu et al. 2011). Around young 
neutron stars, we may see a rotation-powered pulsar wind nebula (Gaensler & Slane 2006). 

In the case of magnetars, the star's persistent X-ray luminosity is much higher than 
its rotational energy loss rate. Since the persistent X-ray luminosity is from magnetic field 
decay, it is possible that a particle fiow (i.e., a magnetism-powered particle wind) is also 
produced during the decay of the star's magnetic field. The luminosity of this particle wincil] 
can be as high as the star's persistent X-ray luminosity, therefore it can also be much higher 
than the star's rotational energy loss rate (Duncan 2000; Section 2.3 below). This particle 
wind will "comb out" the magnetic field lines in the closed field line regions (Harding et al. 
1999). The net result is an enhanced rotational energy loss rate for a given dipole magnetic 
field (Harding et al. 1999; Thompson et al. 2000; Section 3 below). In this "wind aided" spin 
down scenario, the corresponding dipole magnetic field will be much lower than the magnetic 
dipole braking case (Harding et al. 1999; Section 4 below). Wind braking of magnetars will 
also help us to explain recent observations challenging the the traditional magnetar model 
(Section 4 and 5 below). 

Below, we assume that the star's dipole magnetic field is constant during its life time. 
It is the evolution and variation of particle wind luminosity that cause the evolution and 
variation of AXPs/SGRs' timing properties. For example, a short term variation of particle 
wind will cause a variation of the star's period derivative and also contribute to its timing 



^By saying "particle wind", we always mean a mixture of relativistic (or mildly relativistic) particles 
and electromagnetic waves. In Section 3, we will point out the difference of particle luminosity and wind 
luminosity. At present and for general discussions, we will simply use the term "particle wind" . 
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noise. 



2.2. Observational clues for the existence of a particle wind 

The existence of a (rotation-powered) particle wind in normal pulsars is well established. 
The observations of intermittent pulsars give direct support for the existence a particle wind 
(Kramer et al. 2006; Camilo et al. 2012). However, the existence of a magnetism-powered 
particle wind in magnetars is still unknown. Below we will give several observations of AXPs 
and SGRs, which may provide some hints for the existence of a particle wind. 

1. The AXP IE 2259+586 experiences an enhanced period of spin down during outburst 
(Kaspi et al. 2003). Variations of period derivative arc also seen in AXP IE 1048.1- 
5937 (Gavriil & Kaspi 2004), AXP XTE J1810-197 (Camilo et al. 2007), SGR 1806-20 
(Woods et al. 2007) and AXP IE 1547.0-5408 (Camilo et al. 2008) etc. A decreasing 
period derivative is also observed in the radio-loud magnetar, accompanying by a 
decaying X-ray luminosity and radio luminosity (Levin et al. 2012; Anderson et al. 
2012). This may be due to a decaying particle wind during outbursts. 

In the absence of a strong particle wind, an untwisting magnetosphere of a magnetar 
may explain the decreasing period derivative (Beloborodov 2009). However the dipole 
field is of large scale compared with that of multipole field, a varying dipole field 
(especially short time scale variations) is hard to accomplish (Camilo et al. 2007; 
Levin et al. 2012). In the case of wind braking of magnetars, the global dipole field is 
unchanged. Since the particle wind may be the consequences of small amplitude seismic 
activities (Thompson & Duncan 1996), it can vary dramatically even on short time 
scales. A varying particle wind will cause a varying period derivative. The long term 
decay of particle wind luminosity during outburst can account for the decreasing period 
derivative (e.g., AXP XTE J1810-197, Camilo et al. 2007; the radio-loud magnetar. 
Levin et al. 2012). During an outburst, we should expect the wind luminosity first 
increases then decreases. This will cause the period derivative first increases then 
decreases, which may be the case of SGR 1806-20 (Woods et al. 2007) and AXP IE 
1547.0-5408 (before outburst, Camilo et al. 2008). 

2. AXPs and SGRs have a higher level of timing noise than normal pulsars (Gavriil 
& Kaspi 2002; Woods et al. 2002; Archibald et al. 2008). The timing noise may 
be correlated with period derivatives. The timing noise of normal pulsars may be the 
result of a varying (rotation-powered) particle wind (Lyne et al. 2010; Liu et al. 2011). 
Then it is possible that AXPs and SGRs are also braked down by a particle wind. Since 
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AXPs and SGRs are magnetism-powered, the particle wind may also from magnetic 
field decay, i.e., a magnetism-powered particle wind. This magnetism-powered particle 
wind may vary significantly with time, similar to the magnetar's persistent X-ray 
luminosity. Then it may cause a higher level of timing noise in magnetars than that 
in normal pulsars and HBPSRs. 

3. If AXPs and SGRs harbor a strong enough particle wind (either rotation-powered or 
magnetism-powered), then we should see a pulsar wind nebula around the putative 
star. If the the particle wind is magnetism-powered , the same as that of their per- 
sistent X-ray luminosities, then we should see some correlation between the nebula 
luminosity and the stellar luminosity. A possible extended emission is found around 
AXP IE 1547.0-5408 (Vind & Bamba 2009). The luminosity of the extended emission 
is correlated with the star's luminosity (Olausen et al. 2011). Therefore, the extended 
emission around AXP IE 1547.0-5408 may be a magnetism-powered pulsar wind neb- 
ula instead of a duct scattering halo. If this is confirmed in the future, it will be a 
strong evidence for the existence a magnetism-powered particle wind in magnetars. 

A magnetism-powered pulsar wind nebula may also accelerate particles to very high 
energy and radiate high energy photons. An extended emission and a TeV source 
are both seen in the case of SGR Swift J1834.9-0846 (Kargaltsev et al. 2012). If 
the extended emission is found to be a pulsar wind nebula and the associated with 
the TeV source is confirmed, then it is also likely to be a magnetism-powered pulsar 
wind nebula^. A candidate pulsar wind nebula which may contain magnetic energy 
contribution is seen around RRAT J1819-1458 (Rea et al. 2009). 

In summary, there are many uncertainties and ambiguities if we attribute the above 
observations to a particle wind in magnetars. However, we do not know whether AXPs and 
SGRs have a (magnetism-powered) particle wind or not. The possibility of such a particle 
wind can not be ruled out by present observations either. A magnetism-powered particle 
wind in magnetars is helpful to our understanding of the different observational aspects stated 
above. Therefore, the above observational facts may give us some clues for the existence of 
a particle wind in magnetars. Whether a particle wind really exists or not can be tested by 
future studies. 



^After we put this paper on the arXiv (1205.1626), Younes et al. (Jar Xivjl 206. 3330) propose that the 
extended emission of SGR Swift J1834. 9-0846 may be a magnetism-powered pulsar wind nebula (since it 
has a high conversion efficiency). This observation is consistent our analysis here. 
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2.3. Estimation of wind luminosity 

In the magnetar model, the bursts and outbursts are related with the magnetar's seismic 
activities (Thompson & Duncan 1995, 1996). If the observable bursts are associated with 
large amplitude seismic activities, then the low amplitude seismic activities may mainly 
result in a particle wind (Thompson & Duncan 1996). According to Thompson & Duncan 
(1996, eq.(71) there), the particle wind luminosity is 

where is crustal field strength, t is the star's age, and AR^. is the crustal thickness. The 
above equation is only valid for crustal field strength less than 6 x 10^^ G, above which the 
crust may undergo plastic deformations. 

The persistent X-ray luminosity of AXPs and SGRs are from magnetic field decay, e.g., 
internal heating (Thompson & Duncan 1996) or magnetospheric current heating (Thomspon 
et a. 2002; Beloborodov & Thompson 2007). A particle wind may also be produced during 
this process. Since the particle wind and the persistent X-ray luminosity are from the same 
energy reservoir, a natural estimation of the particle wind luminosity is (Duncan 2000) 

Lp~Lx~103^ergs-\ (2) 

which is valid for most AXPs and SGRs. For the transient magnetars, they have a lower 
quiescent X-ray luminosity. Their particle wind luminosity may also be correspondingly 
lower. 

In the wind braking scenario, magnetars are neutron stars with strong multipole fields. 
The strong twisted mangetic field in the vicinity of magnetars will accelerate particles to 
very high energy. Thus, a corona of high energy particle will be formed (Beloborodov & 
Thompson 2007). The footprint of magnetic field lines are anchored to the stellar crust. In 
the presence of frequent low amplitude seismic activities, the corona of magnetars will be 
disturbed continuously. The excitation of such a particle wind in magnetars may be due to 
their seismic activities, especially small amplitude seismic activities (Thompson & Duncan 
1996; Thompson & Duncan 2001; Timokhin et al. 2008). The particles in the magnetar 
magnetosphere can fiow out in two ways. (1) During bursts and giant flares. This burst 
component of particle wind has its duty cycles (Thompson & Duncan 1995; see numerical 
simulations of Parfrey et al. 2012; Section 4.5 below). (2) During persistent state. The long 
term average of many small amplitude seismic activities may result in a persistent particle 
outflow of magnetars (Thompson & Duncan 1996; Duncan 2000). We will mainly focus on 
the persistent component of particle wind. 
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In conclusion, we have already some observational clues for the existence of a particle 
wind in magnetars. Their luminosities can also be estimated, although the underlying mech- 
anism is still lacking. Since both the magnetar's persistent X-ray luminosity and the particle 
wind are from magnetic field decay, the particle wind luminosity may be as high as their 
persistent X-ray luminosities. Therefore, the particle wind luminosity in magnetars can be 
much higher than their rotational energy loss rate. The existence of such a strong particle 
wind will modify the spin-down behavior of magnetars qualitatively. 



The magnetospheres of pulsars and magnetars contain regions of open and closed mag- 
netic field. The closed field lines extend to the hght cyhnder radius in the case of normal 
pulsars (Contopoulos & Spitkovksy 2006). In the case of magnetars, the closed field hue re- 
gion may be smaller. In the presence of a strong particle wind, the natural radial extension of 
closed field line regions is the radius where the kinetic energy density of particle wind equals 
the magnetic energy density (Harding et al. 1999; Thompson et al. 2000). Particle flows in 
the closed field line regions belong to the domain of closed field line region electrodynamics 
of magnetars (Thomspon et al. 2002; Beloborodov & Thompson 2007; Tong et al. 2010b). 
Particle flow coUimated around the polar cap may dominate the spin down of the central 
star. The opening angle of the polar cap region is determined by the coupling between the 
magnetar crust and its magnetosphere. The total particle luminosity Lp is determined by 
the decay of magnetic field energy. Only a fraction of this particle wind can flow out to 
infinity and contribute to the spin down of the magnetar. The escaping particle luminosity 
is denoted as L„, i.e., wind luminosity. Then it is natural that L„ < Lp. For a given particle 
luminosity, the maximum braking case is accomplished when the wind luminosity equals the 
total particle luminosity. 



For a neutron star with angular velocity fl = 27r/P (P is rotation period), its light 
cylinder radius Ric is (the radius where the rotational velocity equals the speed of light) 



3. Rotational energy loss rate due to a peirticle wind 



3.1. Description of the global magnetospheric structure 



3.2. The simplest case: — Lp 




(3) 
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where c is the speed of hght. In the case of magnetars, with the aid of a particle wind, 
the magnetic field hnes are combed out at a radius Topen (where the particle energy density 
equals the magnetic energy density, Harding et al. 1999) 

/ r2„2„\ 1/4 / r2^2„\ 1/4 

ro„e„ = r„(^) =r„(^) = 4.1 x 10»(,;%l^\m, (4) 

where ro = 10^ cm is neutron star radius, Bq = bo x Bqeb is dipole magnetic field at the 
magnetic pole, and = Lp = i^p,35 x 10^^ erg s^^ is the particle wind luminosity (assuminjf] 
L„ = Lp, and assuming the escaping particle wind becomes near isotropic at ropcn)- The 
polar cap radius now is 

Rpc = ro(ro/ropen)'/' = 1.6 x cm. (5) 

The corresponding polar cap opening angle is 

Cn = ^o/ropen = 2.4 X 10"^ 6"^/' (6) 

Typically, 6'open = 1-6 x 10^^ 6q ^^"^ Lp^^^. The polar cap opening angle 6'open depends on the 
wind luminosity L^. 

This forms the basic structure of a wind-loaded magnetosphere. The star may form 
a current circuit in the open field line regions. The rotational energy loss rate due to this 
particle wind is (Harding et al. 1999) 



'"open 



2 



For the traditional magnetic dipole braking, the corresponding rotational energy loss rate 
i£| (Shaprio & Teukolsky 1983) 

E, = ^fi^. (8) 

Therefore, eq.(I7]) can be rewritten as 

2 ^ /L/--'" 



■^This means that there is only a small particle flow in the closed field line regions. 

Note that eq.® is for orthogonal rotators. While in the wind braking case, for simplicity, we are 
considering an aligned rotator (Harding et al. 1999). Therefore, eq. ^ should be taken as definition rather 
than derivation. 
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A second way to calculate the rotational energy loss rate due to a particle wind is 
provided by Thompson et al. (2000)0. The outflowing particles will corotate with the star 
up to the radius ropen- For relativistic (also mildly relativistic) particles, the rotational energy 
carried away by this particle wind is (Thompson et al. 2000) 

^w = ?VroU = 4^dfi^V''. (10) 



3c2 ^ "Vi^d/ 

A third way to calculate the rotational energy loss rate due to a particle wind can be 
done in analogy with that of Xu & Qiao (2001). The electric current in the two polar caps 
will carry away the rotational energy of the star in the presence of an acceleration potential. 
This acceleration potential is due to unipolar induction. Assuming maximum acceleration 
potential, the rotational energy loss rate is 

1 /2 

K = 2/pc$max = ^^d (^J-^ , (11) 

where Ipc = nRp^.pGjc is the polar cap current (for one polar cap), pcj is the Goldreich- 
Julian density, and $max is the maximum acceleration potential due to unipolar induction 
(Ruderman & Sutherland 1975) 

*^-^(^)^- 

Therefore, irrespective of the details of particle wind, accurate within a factor of two, 
the rotational energy loss rate due to a particle wind can be written as (Harding et al. 1999) 

r \ 1/2 / r \ 1/2 



E. = E,{fJ =E.[fJ (13) 



The second identity is obtained by assuming = Lp. 
From equation f|T3|) . we see that 



1. For a rotation-powered particle wind Lp ~ —E^ot, ~ -Ed ~ —E^ot, wind braking is 
quantitatively similar to the case of magnetic dipole braking in vacuum. The effects 
of particle wind will mainly cause higher order modifications, e.g., a different braking 
index etc. This is the case of normal pulsars. 



^''opcn in Harding et al. (1999) is equivalent to the Alfven radius i?A in Thompson et al. (2000), except 
for a difference of constant factor 2^/^. Hereafter, ropon is employed. 
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2. For magnetars, there may be a magnetism-powered particle wind Lp ^ — -Erot- Wind 
braking of magnetars will result in = — -Erot ^ E^. Therefore, the magnetic 
dipole braking is enhanced due to the presence of a magnetism-powered particle wind 
(Harding et al. 1999). This will cause a strong reduction of magnetar's dipole field. 
Meanwhile, high order effects will also exists, e.g., a different braking index, larger 
timing noise, a magnetism-powered pulsar wind nebula, etc. 



3.3. Detailed considerations of wind luminosity 

The above simplest case assumes the escaping wind luminosity is equal to the total 
particle luminosity. From equation (|6]), the polar opening angle depends on the escaping wind 
luminosity. This means that the polar cap opening angle (at the star surface) is affected by 
the physics happening at ropen- It is not known how this is accomplished. Alternatively, the 
polar cap opening angle of the particle wind may be an independent parameter. The total 
particle luminosity may involve a particular angular distribution. This angular distribution 
may result from coupling between the magnetar crust and its magnetosphere. The typical 
time scale of this coupling may be estimated from quasi-periodic oscillations in magnetars 
(Timokhin et al. 2008; Watts 2011). The fundamental frequency is about u ~ 20 Hz. The 
length scale of coupling between the neutron star and its magnetosphere is 

c o/20Hz\ , , 

rmax ~ — ~ 5 X loM cm. 14 

Su \ u J 

The corresponding polar cap opening angle is 

=^~2xl0-=^f^). (15) 



.20 Hz. 

1 /2 

Typically, 6^ ~ 0.05 {2^^) ■ The particles will mainly flow through the polar cap area 
with opening angle 6s . In the following calculations, we will take 6s as the fundamental 
input parameter, rmax etc will be functions of 6s. 

The particles from the two polar cap regions can flow out to radius larger than rmax- 
Considering the presence of strong magnetic field, a significant amount of the outflowing 
particles may be trapped in the closed field line regions in the magnetospher^. Only a 
fraction of them can flow out to infinity and therefore carry away the star's rotational energy. 
The Alfven radius characterize the effect of magnetic field quantitatively. We denote it as 



^These trapped particles may contribute to the persistent X-ray emissions of magnetars. 



- 12 - 



Topen in accordance with equation (jl]). In the present case, it is also defined as the radius 
when the particle energy density equals the magnetic energy density 

ip{r)c' ~ (16) 

where 7 and p(r) are the Lorentz factor and mass density, respectively. When particles move 
along magnetic field lines, their kinetic energy is conserved (not considering radiation losses). 
The mass density may scale with the local Goldreich- Julian charge density p(r) oc pcj oc 
Therefore 

where ps is the mass density at the star surface. According to the definition of particle 
luminosity and assuming uniform distribution across the polar cap region 

Lp = 27r(ro0s)'7Psc'c, (18) 

then Topen is 



ropen - '^O^ p (iyj 



Svr Lp 

= 7x lO^foo/'Lp Jf (0s/O.O5)'/' cm. (20) 

Only the escaping wind particles can carry away the star's rotational energy. From 
the definition of wind luminosity, oc 0"^^^^ oc 1/ropen- At the same time, the total par- 
ticle luminosity is (x 6l (x 1/rmax- The wind luminosity is related to the total particle 
luminosity 

= l/-^. (21) 

^open 

Taken the polar cap opening angle as the fundamental parameter, rmax will be rmax = fo/Gl- 
Therefore, the wind luminosity is 

L„ = 6 X 10^3 6-'^'^pS(^s/0.05)-«/3 erg s-\ (22) 

The wind luminosity depends strongly on the polar cap opening angle, i.e., how the neutron 
star couples with the magnetosphere. In the present case, the wind luminosity is a fraction 
of the total particle luminosity. Then, it must be that < L^. In terms of r^ax and ropen, 
it must be that rmax < ''^open- 

The calculation of rotational energy loss rate is the same as the previous section. From 
equation (|T3l) . the rotational energy loss rate due to a particle wind in the present case is 

/r \i/2 

E^ = eJ^] , (23) 
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where is determined by equation (l22l) . The neutron star's dipole magnetic field is obtained 
by equahng -E,ot = E^, 

Bo = 3.3 X 10 
= 3.3 X 10 

The dipole magnetic field is determined by four parameters: the period and its derivative, 
the total particle luminosity, and the polar cap opening angle. If the polar cap opening angle 
is three times smaller, the dipole magnetic field will be ten times lower. 

In conclusion, considering detailed modeling of wind luminosity, the rotational energy 
loss rate is reduced compared with the simplest case. The model parameter space is larger 
with the addition of another variable 6s- There are parameter space that the corresponding 
dipole magnetic field is only slightly lower than the magnetic dipole braking case. At the same 
time, there are also some parameter space that the dipole magnetic field is much lower than 
the magnetic dipole braking case. The following calculations in Section 4 are mainly done 
in the simplest case. This corresponds to maximum braking for a given particle luminosity. 
In this way, we want to demonstrate to which extent can wind braking of magnetars help 
to explain the current observations. For the calculations in Section 4.2, 4.4, and 4.5, the 
conclusions are unaffected by different assumptions. For the calculations in Section 4.1 and 
4.3, the results may only change quantitatively. 



P 



3/2 



75 L-Ue,/o.o5fG 




(24) 



4. Wind braking of magnetars 

Wind braking of magnetars had been considered previously by Marsden et al. (1999, for 
the case of SGR 1900+14), Harding et al. (1999, for the case of SGR 1806-20), Thompson 
et al. (2000, for the case of SGR 1900+14). They mainly talked about wind braking during 
outbursts, although some of the formulae of the long term wind-aided spin down are also 
given by Thompson et al. (2000, Section 4.1 there). We explore the wind braking in more 
details and apply it to all magnetars. A comparison with up-to-date observations is also 
presented. 
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4.1. Dipole magnetic field 

From eq. (j4]) and eq. ( 1T3|1 . the presence of particle wind amplifies the magnetic dipole 
braking rotational energy loss rate. Therefore, wind braking is valid only when 

Lp>Ed. (25) 

Equaling the rotational energy loss rate —Ej-ot (= —I^^) and eq. f|T3l) . we get 

- Erot = < L^- (26) 

Wind braking of magnetars is valid only when the wind luminosity is greater than the star's 
rotational energy loss rate. Equation (|26|) can be rewritten as 

— -E'rot = E„ > E^. (27) 

The characteristic magnetic field obtained by assuming magnetic dipole braking is only the 
upper limit of the star's true dipole magnetic field. 

Assuming magnetic dipole braking 

- -E'rot = -E'd = — ^ — , (28) 
the dipole magnetic field (at the magnetic pole) is 

Bo = 6.4 X lO^^/p^G = 6.4 X lO^M ^1 G. (29) 

" \^10slO-ii J ^ ' 

It is two times larger than usually reported since the polar magnetic field is two times larger 
than the equatorial magnetic field (eq. (5.17) in Lyne & Graham-Smith 2012 and corre- 
sponding discussions). However, the above magnetic dipole braking is originally designed 
for rotation-powered pulsars. Magnetars may be wind braking instead of magnetic dipole 
braking, as discussed above. In the case of wind braking 

/L 

- 4ot = = . (30) 



E, 



The corresponding dipole magnetic field is 

So = 4.0 X 10^^^ L^lf G = 4.0 X 10^^^^ L'J^ G. (31) 



For typical AXPs and SGRs, the dipole magnetic field in the case of wind braking is about 
ten times lower than that of magnetic dipole braking. Therefore, AXPs and SGRs may 
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be magnetars without strong dipole field. Only strong multipole field (~ 10^^ — 10^^ G) is 
required to power their bursts, persistent emissions, and braking. 

At the time when Harding et al. wrote their wind braking paper (Harding et al. 1999), 
they did not realize that there are two kinds of magnetic fields in magnetars: dipole field 
and multipole field. When they saw that a strong dipole field is not needed in the case of 
winding braking, Harding et al. said that "the magnetar model must be abandoned" as 
the penalty of wind braking. With the presence of multipole field, AXPs and SGRs can 
also show magnetar-like activites without a strong dipole field. This point is demonstrated 
clearly by the observation of SGR 0418+5729 (Rea et al. 2010). The timing of SGR Swift 
J1822.3— 1606 further strengthens this point (Rea et al. 2012a). 

Table [T] summaries the observed parameters and deduced quantities for all AXPs and 
SGRs (17 in total), which have period, period derivative, and persistent X-ray luminos- 
ity measured. Figure [1] shows the magnetar persistent X-ray luminosity versus the star's 
rotational energy loss rate. We employ the following two ways to the model the particle 
luminosity from magnetars. 

1. Considering that for all AXPs and SGRs, they must have a strong multipole field 
(~ 10^^ — 10^^ G) in order to show magnetar-like activities. This is also true for low 
magnetic field SGRs (Rea et al. 2010, 2012a). Therefore, if the total field strength 
determines the particle luminosity, the particle luminosity will be more or less the same 
for all magnetars. In this case, we assume a particle luminosity Lp = 10^^ erg s~^ for 
all sources. From Figure [T|, we see that all AXPs and SGRs are braked by a particle 
wind except AXP IE 1547.0-5408. For AXP IE 1547.0-5408, the effect of a particle 
wind will mainly result in high order spin down behaviors, e.g. a magnetism-powered 
particle wind surrounding the putative magnetai0. 

2. On the other hand, different sources may have a different evolution history. Irrespective 
of the detailed wind mechanism, the magnetar's particle luminosities may follow their 
persistent X-ray luminosities. In this case, we assume that the particle luminosities 
are the same as their persistent X-ray luminosities. From Figure [H except for the five 
sources with < — -Erot, the rest of AXPs and SGRs are all braked down by a particle 
wincPI. 



^For the case of AXP IE 1547.0-5408, since it rotational energy loss rate is also relatively large, the 
surrounding pulsar wind nebula may be a mixture of rotation-powered and magnetism-powered particle 
wind. 

^this may explain the "fundamental plane" of magnetar radio emission, see Section 4.2 below. 
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At present, we do not know the detailed mechanism of magnetar wind. The actual case may 
lie between these two extremes. 

Figure [2] and [3] show the dipole magnetic field in the case of wind braking versus the 
dipole magnetic field in the case of magnetic dipole braking. From Figure [2] and [3l we see 
that 

1. For most AXPs and SGRs, their dipole magnetic field by assuming wind braking are 
ten time lower than that of magnetic dipole braking. This may help us to understand 
why the magnetar supernova energies are of canonical value (Vink & Kuiper 2006; 
DairOsso et al. 2009). 

Numerical simulation of particle wind during magnetar bursts also suggests that the 
long term averaged period derivative may be greatly amplified (Parfrey et al. 2012). 
The actual dipole magnetic field may be significant lower than the magnetic dipole 
braking case. This is consistent with our considerations here. 

2. The corresponding dipole magnetic field -Bdip.w ranges from 10^^ G to 10^^ G. A strong 
dipole magnetic field (Kouveliotou et al. 1998, > -Bqed = 4.4 x 10^'^ G) is no longer 
a necessary input. In the wind braking scenario, magnetars are neutrons with strong 
multipole field. For most sources, their dipole field may or may not be as strong as 
their multipole field. 

3. For several sources, their -Bdip,w are in the range 10^'^ G — 10^'^ G. This is similar to 
that that of X-ray dim isolated neutron stars (Kaplan & van Kerkwijk 2011; Tong et 
al. 2010b). Therefore, when the magnetar activities of these sources calm down, they 
will become X-ray dim isolated neutron stars naturally. 

4. There are now more low magnetic field magnetars, with -Bdip,w < 4.4 x 10^^ G. There- 
fore, in the case of wind braking, SGR 0418+5729 (Rea et al. 2010) is not such peculiar 
as before. 

Furthermore, from eq. (IHUj) and eq. (I3T1) we see that for a given magnetar 

1. The variation of particle wind will result in a variation of P, P oc Lp^^. The may explain 
the variation of P of many AXPs and SGRs (see Section 2.2 above). A decreasing 
particle wind will result in a decreasing period derivative during magnetar outbursts 
(Camilo et al. 2007; Levin et al. 2012; Anderson et al. 2012). 

2. Although magnetars and high magnetic field pulsars (HBPSRs) are close to each other 
on the P — P diagram, they may be totally different from each other. In the case of 
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wind braking, magnetars are neutron stars with strong multipole field. While HBPSRs 
may be neutron stars only with strong dipole field. 

Observationally, AXPs and SGRs have a larger level of timing noise (Gavriil & Kaspi 
et al. 2002; Woods et al. 2002; Archibald et al. 2008). This may be the result that 
they are wind braking instead of magnetic dipole braking. Meanwhile, most HBPSRs 
do not show magnetar-like activities which may be that most of them do not have as 
strong multipole fields as magnetars (Ng & Kaspi 2011; Pons & Perna 2011)1^ 

The above calculations are done by assuming L„ = Lp. As discussed in Section 13. 3[ 
a strong reduction of dipole magnetic field is possible only when L„ is comparable to Lp. 
This corresponds to a very coUimated particle wind at the star surface. A small reduction 
of dipole magnetic field results from detailed modeling of wind luminosity. Assuming a 
constant polar cap opening angle = 0.05 and Lp = L^, the corresponding dipole magnetic 
field is also shown in tableland Figure H] (The case is similar when assuming = 0.05 and 
Lp = lO^^ergs"^). As pointed out in the beginning of this section, if the wind luminosity 
is lower than the rotational energy loss rate, the dipole magnetic field will be the same as 
that of magnetic dipole braking. Then, the effect of a particle wind will mainly reflected in 
higher order modiflcations, e.g. braking index, timing noise, a magnetism-powered pulsar 
wind nubulae etc. Considering detailed modeling of wind luminosity, the dipole magnetic 
fleld will be the same as the magnetic dipole braking case for most sources. Only for four 
sources, their dipole magnetic flelds are lower than the magnetic dipole braking case. 



4.2. Acceleration potential 

Most of the electromagnetic emission of magnetars is thought to originate in the closed 
fleld line region (Thompson et al. 2002; Beloborodov & Thompson 2007; Tong et al. 2010b). 
Meanwhile, since the rotational energy is always present, we should also see some rotation- 
powered activities in magnetars (Zhang 2003). The rotation-powered activities are almost 
inevitable especially when we assume that AXPs and SGRs are also magnetic dipole braking 
as rotation-powered pulsars (Tong et al. 2011). The acceleration potential in open fleld line 
regions characterizes this point quantitatively. 

The maximum acceleration potential in pulsar open fleld line regions is (Ruderman & 



^If a HBPSR also has a strong multipole field, it can also show magnetar-like activities. This may be the 
case of PSR J1846-0258 (GraviU et al. 2008; Pons k Perna 2011). 



- 18 - 



Table 1: Measured quantities and inferred dipole magnetic field of magnetars. 
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PSR J1622-4950 


4.33 
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0.0063 
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-Sdip^d 
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Notes: Column one to eight are respectively, source name, period, period derivative, persistent X-ray 
luminosity in the 2 — 10 keV range, dipole magnetic field assuming magnetic dipole braking, dipole 
magnetic field in the case of wind braking assuming a wind luminosity = Lp = lO^^ergs"^, dipole 
magnetic field in the case of wind braking assuming a wind luminosity = Lp = L^, and dipole magnetic 
field in the case of wind braking assuming 9s = 0.05 and Lp = L^- All data are from the McGill SGR/AXP 
online catalogue (http://www.physics.mcgill.ca/~pulsar/magnetar/main.html, up to January, 27, 2012), 
except for SGR 0418+5729 (from^Rea et al. 2010), and Swift J1822.3-1606 (from Rea et al. 2012a). The 
first column are ordered roughly by the source's discover time. 

median value is used during calculations 
'': upper limit is used during calculations 

*: When the wind luminosity is smaller than the rotational energy loss rate, the dipole magnetic field in 
the case of wind braking is the same as that of magnetic dipole braking. A magnetism-powered particle 
wind mainly results in a magnetism-powered pulsar wind nebula and other high order modifications. See 
text for details. 



- 19 - 



10" 



10" 



10" 



10 








■ 

4U0142 

■ 

1E2259 


1 1 

1E1841 ■ ■<;rpisi(fe 

J1708J0100 / 
/ ■ 

/ J1714 








■ / 


■ 

SGR1627 




■ 

SwiftJ1822 


■ 

J1647 / 


/ ■ 

/ J1622 


■ 

1E1547 


■ 

SGR0418 






■ 

J1810 

1 





29 


10^' 




10" 


10^'^ 



-£iot(erg s ) 



Fig. 1. — Persistent X-ray luminosities of magnetars versus their spin down luminosities. 
The solid line is Lx = —Emt- See table [Hand text for details. 



Sutherland 1975) 



$_ = ^f^) . (32) 

In the case of magnetic dipole braking R^^ = ro(ro/-Ric)^^^, the corresponding acceleration 
potential is 

3 —Ey-ot 



1/2 



(33) 



In the case of wind braking, the polar cap radius is given by eq.(l5]). Although the polar cap 
radius is larger than the magnetic dipole braking case, the dipole magnetic field is lower. 
The net effect will be concealed. The corresponding acceleration potential is 



1/2 



rot 



(34) 



The maximum acceleration potential is the same (within a factor of two) in the wind braking 
case and magnetic dipole braking case. 

Although the maximum acceleration potential is the same, the detailed acceleration 
mechanism will be qualitatively different. In the presence of a particle wind, vacuum gaps 
may not be formed, e.g., outer gap etc. This may explain the conflicts between outer gap 
model in the case of magnetars and Fermi observations (Tong et al. 2010a, 2011). Meanwhile, 
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Fig. 2. — Dipole magnetic field in the case of wind braking versus dipole magnetic field in 
the case of magnetic dipole braking. A wind luminosity Lp = 10^^ erg s~^ is assumed for all 
sources. The solid, dashed, and dotted lines are for -Bdip.w = -Sdip,d, O.lSdip.d, 0.01-Bdip,d, 
respectively. The dot-dashed line marks the position of quantum critical magnetic field 
Bqed = 4.4 X 1013 G. 
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Fig. 3. — The same as Figure |2l The wind luminosities are assumed to be the same as their 
persistent X-ray luminosities. See text for details. 



- 22 - 



100 



o 



0.01 





, 1 , , , , 1 , , , 1 , , , , 


, 1 , , 

^/SGR1806 






I3II9OO ,-7 

^5^41 






y 




^/J1647 






/ ^4-00142 




: ■ 


1522S9 




SGR0418 














. 1 , ... 1 , , , 1 . , , , 


. 1 , , , . 



0.5 1.0 5.0 10.0 50.0 100.0 

Sdip.^/CIO^G) 
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space charge limited flow type acceleration mechanism may still exist (Xu 2007). In a wind 
loaded magnetosphere, detailed calculations of space charge limited flows are needed in the 
future. 

In calculating Figure [3l we show that only those sources with > —E^ot are wind 
braked down. While for sources with < — -Erot, they are still magnetic dipole braking, 
the same as rotation-powered pulsars. A magnetosphere similar to that of rotation-powered 
pulsars is prepared during the persistent state. This may be taken as the initial state. An 
outburst will may trigger the radio emission of magnetars as observed. Then it is natural that 
only sources with < — -E'rot can have radio emissions. This may explain the "fundamental 
plane" of magnetar radio emissions found by Rea et al. (2012b). More detailed investigations 
are needed. 



4.3. Spin down evolution and age 

A given magnetar, with dipole magnetic field Bq = bo x -Bqed = bo x 4.4 x 10^^ G 
{bo ~ 1 from eq. fl3Tl) ). and a wind luminosity = Lp = i^p,35 x lO^^ergs"^, will evolve 
from magnetic dipole braking at early stage to wind braking at later stage. At present, we 
assume that Bo and Lp are both constants (A decaying particle wind will be considered in 
Section 5.1 below). From eq. fl2^ and eq.®, at the early stage, the star rotates very fast and 
Ed is larger than Lp. Therefore, the star will be braked down by magnetic dipole radiation 
at early stage. However, at later stage, the star will be slowed down and Ed will be smaller 
than Lp. Therefore, the star will become wind braking at later stage. The transition from 
magnetic dipole braking to wind braking happens at 

= = (35) 

The corresponding rotation period is 

Pi=0.66 6fL-j//s. (36) 

Pi can also be obtained by requiring ropen < -Ric (Thompson et al. 2000). When the star's 
rotation period is less than Pi it will be braked down by magnetic dipole radiation. The 
corresponding period derivative at the transition point is 

P, = 7.2xW~'HfLl%. (37) 

If the magnetar rotation period at birth is much less than Pi, then the star age at Pi is 

h = Tci ^ ^ = 1.4 X lOHo-iLpjf yr. (38) 
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The transition age ti is similar to the supernova remnant age associated with AXPs 
(Vink & Kuiper 2006). Beginning from ti, Lp > — -Erot (eq.f l26p ). the star will be braked 
down by a particle wind. Furthermore, the particle wind of magnetars is from magnetic 
energy decay Lp ~ —E-^, where is the star's magnetic energy stored mainly in the 
form of multipole field. Therefore, during wind braking phase, —Eb > — -Erot- The star's 
activities will be dominated by magnetic energy output rather than rotational energy output. 
AXP/SGR-like activities may appear, i.e., the pulsar becomes a magnetar. 

We now consider how a magnetar evolves from {Pi, Pi) to {P2,P2) (Thompson et al. 
2000). When we assume Bq and Lp are both constants, then from eg. (1301) . at wind braking 
phase 

P Pi P2 

— = — = — = constant. (39) 
P Pi P2 

The period will evolve with time as 

P2 = Piexp{^^}, (40) 
2rc,i 

where ^2 and ti are the star's true age at P2 and Pi, respectively. Tc^i is the characteristic 
age at Pi. For transition from magnetic dipole braking to wind braking, ti = r^^i. However, 
in the general case, ti is not always equal to Tc^i. The star's age at a given period Pj is 

t2 = ti + 2r,,ilog^. (41) 

After ti, the star's period increases exponentially. For P2 not very large than Pi, we have 
t2 ~ ti = Tc^i = Tc^2, where Tc,2 is the star's characteristic age at P2. 



4.4. Braking index 

The braking index of a pulsar is defined as (Shapiro & Teukolsky 1983) 

n = - (constant )^]", (42) 

where n is called the braking index, n = 3 for magnetic dipole braking. For wind braking, 
from eq. (15CT]) , we have 



6r)4\ 1/2 



Therefore, n = 1 for wind braking (assuming Bq and Lp are both constants). The braking 
index of PSR J1734-3333, n = 0.9 ± 0.2, may imply that it is of wind braking (Espinoza 
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et al. 2011; a rotation-powered particle wind). Future braking index measurement of a 
magnetar will help us make clear whether magnetars are magnetic dipole braking or wind 
braking. Because the braking index will deviate from one if Bq and/or Lp changes with time, 
a braking index of a magnetar may also tell us the evolution of its particle wind. 

4.5. Duty cycles of particle wind 

Harding et al. (1999) considered the duty cycles of a particle wind whose luminosity 
is Lp ~ 10^''ergs~^. It is shown that, due to the duty cycles of particle wind, the dipole 
magnetic field and age vary continuously from the dipole braking case to the wind braking 
case (Figure 1 in Harding et al. 1999). However, the particle luminosity considered by 
Harding et al. (1999) is much stronger than we considered here Lp ~ lO^^ergs"^. It is 
possible that there are two types of particle wind: 

1. A persistent component associated with the magnetar 's persistent emissions. The 
particle luminosity is Lpp ~ Lx ~ lO'^^ergs"^. 

2. A burst component associated with outbursts of magnetars. The corresponding particle 
luminosity may be about Lpb ~ Lburst ~ lO^^ergs"^. 

The burst component of particle wind may contribute to the enhanced spindown of magnetars 
after glitches (Kaspi et al. 2003) and the possible "radiation braking" during giant flares of 
SGR 1900+14 (Thompson et al. 2000; Parfrey et al. 2012). 

The long term averaged spindown of magnetars can be modeled similarly to that of 
Harding et al. (1999) 

~ {Erot) = -£'w, burst -Dp + -E^, persistent ( 1 ~ -^p)? (44) 

where Dp is the duty cycle of the burst component of particle wind. From eq. fll3p . the above 
equation can be rewritten as 

- (4ot) = E'/'LlCDp + eTlI\1 - Dp) = eH'l]!^ (45) 

1 /2 1 /2 1 /2 

where Lj^ = L^^^ Dp + Lpp {1 — Dp) is the effective particle luminosity. For typical param- 
eters, the effective particle luminosity is 

= 10 L%7D, + Ll!p%{l - Dp). (46) 

For Dp = 0, this is just the case we considered above. For Dp = 1, this corresponds to a 
strong wind case (Lp = 10^''ergs~^) as considered by Harding et al. (1999). The duty cycle 
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can be estimated from the observations of the transient magnetar SGR 1627-41 (Mereghetti 
et al. 2009). The duration between two outbursts is about ten years. Therefore, the 
maximum value of duty cycle is about 0.1. The corresponding effective particle luminosity 
is lIq^^ = -Z^pb^37 + 0.9Lpp^35, about two times larger than the persistent component of the 
particle wind. In conclusion, the previous discussions are still valid considering the possible 
existence of a burst component of particle wind. 



5. Discussions 

5.1. A decaying particle wind 

In the magnetar model, both the persistent and burst emissions of AXPs and SGRs are 
powered by magnetic field decay. The total magnetic field will decay with time. Meanwhile, 
the photon luminosity as well as the particle luminosity will also evolve with time. Eventually 
both the photon luminosity and particle luminosity will also decay with time (TuroUa et al. 
2011). In the case of a decaying particle wind, the spin down evolution of magnetars will 
be different from previous considerations. Considering different avenues for magnetic field 
decay, the total magnetic field may decay with time in a power law form (Heyl Sz Kulkarni 
1998). The consequent magnetic energy decay rate —Eb will also of power law form. Since 
the particle luminosity is from the magnetic energy decay, we may assume a power law form 
of particle luminosity 

L^t) = Lp,o (J-^ , < a < 2, (47) 

where Lp^o and a are constants, to is the time when the magnetic field starts to decay 
significantly, to may be of the same order as ti when wind braking starts to operate. For a 
larger than two, Lp{t) decays more rapidly than E^. In this case, the wind braking criterion 
is not fulfilled (eq. (125|1 ). In the case of decaying particle wind, by integrating eq. (l30!l . we 
can get the spin down evolution of magnetars. For < a < 2, the period evolves with time 

P2 = Pi exp{ ]' '-}, (48) 

and the star age at a given period is 



t2\^fj =ti + (2-a)r,,ilog-^. (49) 
For the special case of a = 2, the corresponding expressions for period and age are 

+ \ ii/2Tc. 



^2 = Pi ( ^ ) , (50) 
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7-) \ 2Tc,l/tl 

-r 2 



h = h [J-J . (51) 

Equation ( 150|) is now the same as the magnetic dipole braking case by setting ti = Tc^i (eq. 
(5.18) in Lyne & Graham-Smith 2012). 



5.1.1. Calculation of braking index 

The braking index predicted for the most luminous AXP (4U 0412+61, Dib et aL 2007) 
is shown as function of age in Figure O -^^p^o = 10^''ergs~^ is assumed. For a constant 
particle wind, the braking index n = 1 is obtained, as previously discussed. For the critical 
case a = 2, the braking index n = 3 is obtained the same as the magnetic dipole braking 
case, as can be seen from eq. (150!) . For the intermediate case < a < 2, a braking index 
n = 1 — 3 is obtained. Future braking index measurement of this source may tell us whether 
it is wind braking or magnetic dipole braking. 
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Fig. 5. — Braking index in the case of wind braking as a function of age. The parameters 
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lines are for a = 0, 0.5, 1, 1.5, 2, respectively. 



5.2. The presence of a fallback disk 

In the case of wind braking, the star's true age is of the same order of the character- 
ristic age t ~ t^. For those magnetars whose supernova remnant age tsnr ~ Tc, then it is 
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understandable that they are wind braking. However, for AXP IE 2259+586, its supernova 
remnant age tsnr ~ 10^ yr ^ Tc = 23 x 10^ yr (Vink & Kuiper 2006). For a decaying par- 
ticle wind, the star's true age can be less than Tc- However, Lp(tsnr) will be larger than 
Lx ~ lO^^ergs"^. Therefore, additional torque may be needed for AXP IE 2259+586. 

The presence of a fallback disk may help to solve this age discrepancy (Shi & Xu 2003). 
At early phase, a fallback disk provides the braking torque of the magnetar. At the end of 
disk braking, the star has been slowed down significantly, e.g. ti = 2 x 10^ yr. Pi = 6.7 s. For 
a particle luminosity Lp = lO'^^ergs"^, the evolution of rotation period is shown in Figure 

El 

Observationally, there may be a debris disk around IE 2259+586 (Kaplan et al. 2009). 
If we assume that SGR 0418+5729 is also a young magnetar, then a fallback disk is also 
needed (in the early stage) to spin down it to the present period (Alpar et al. 2011). For 
the disk torque to operate effectively, the dipole magnetic field can not be too high, e.g., 
Sdip = 10^2 - 10^=^ G is required (Shi & Xu 2003; Alpar et al. 2011). This is consistent with 
the dipole magnetic field obtained by assuming wind braking (see eq. ipT]) ). 

In conclusion, there may be a fallback disk in the early stage of a magnetar. This 
fallback disk may help to solve the age discrepancy. At present, they have been slowed down 
significantly and have become wind braking. 




2000 5000 1x10* 2x10* 5x10* 1x10= 

age (yr) 



Fig. 6. — Evolution of rotation period as a function of age, calculations for AXP IE 
2259+586. The star is AXP IE 2259+586, tsnr is taken as the true age. 
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5.3. Spin down evolution of newly born magneteirs 

Magnetars are thought to be descendants of rapidly rotating proto-neutron stars, with 
rotation period ~ 1ms (Duncan & Thompson 1992). A strong dipole field (-Bdip ~ 10^^ G) 
will cause the spin down time scale of the magnetar less than the supernova shock breakout 
time. This will cause the supernova associated with magnetar birth more energetic (Duncan 
& Thompson 1992). However, studies of supernova remnants associated with AXPs and 
SGRs show that the putative supernova energies are of canonical value (Vink & Kuiper 
2006). This provides challenges to the traditional magnetar model. If magnetars are wind 
braking instead of magnetic dipole braking, then they will have much weaker dipole field. 
The corresponding spin down time scale will be much longer than the shock breakout time, 
Tsd ~ 605^^ (Pi/1 ms)^ hr. This may explain the observations of Vink & Kuiper (2006). 

Moreover, in the presence of strong multipole field, magnetars are of prolate shape. This 
may cause them to emit strong gravitational waves after birth (Dall'Osso et al. 2009). The 
gravitational waves will also carry away some amount of the initial rotational energy. For 
gravitational wave to operate effectively, its competing process (i.e., magnetic dipole braking) 
can not be too strong. Therefore, a weaker magnetic dipole field is required i?dip ^ 10^"^ G. 
This is also consistent with the result of wind braking. In the actual case, a combination of 
these two processes, i.e., longer spin down time scale and gravitational wave emissions, may 
account for the observations. Their contributions depends on the dipole and multipole field 
strength of the star, which may vary from source to source. 

5.4. Magnetism-powered pulsar wind nebula 

A particle wind with luminosity lO^^ergs"^ may produce a visible nebula around the 
central magnetar (The putative nebula may also contain contributions from a rotation- 
powered particle wind). This pulsar wind nebula is magnetism-powered in nature since 
the particle wind is originated from magnetic field decay. There may be a pulsar wind 
nebula around AXP IE 1547.0-5408 (Vink & Bamba 2009). Since both the particle wind 
and the persistent X-ray luminosity of a magnetar are from magnetic field decay, there will 
be a strong correlation between them. Then for a magnetism-powered pulsar wind nebula, 
we should see a correlation between the nebula luminosity and the magnetar luminosity. 
This is just the case of Figure 2 in Olausen et al. (2011). In Olausen et al. (2011), they 
see a strong correlation between the extended emission of AXP IE 1547.0-5408 and its 
source flux. Therefore, Olausen et al. concluded that the pulsar wind nebula origin for the 
extended emission is ruled out and it is a dust scattering halo. However, a strong correlation 
between the extended emission and the source flux just rules out the rotation-powered pulsar 
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wind nebula hypothesis. Such a correlation is a natural result if the pulsar wind nebula is 
magnetism-powered. Future multiband observations of this source may tell us whether it is 
a magnetism-powered pulsar wind nebula or a dust scattering halo. 

For a magnetism-powered pulsar wind nebula, an extreme case is that the nebula lu- 
minosity can exceed that of the star's rotational energy loss rate Lp„n > — -f'rot- However, 
for young magnetars, their rotational energy loss rates are also very high. Therefore, the 
extreme case may be very hard to achieve. A possible case is that we can see a high conver- 
sion efficiency of the putative nebula. The possible pulsar wind nebula seen around RRAT 
J1819— 1458 has a relatively high conversion efficiency (Rea et al. 2009). It may contain 
contributions from a magnetism-powered particle wind. 



6. Conclusions 

Considering recent observations challenging the traditional magnetar model (neutron 
stars with both strong dipole field and strong multipole field), we explore the wind braking 
of magnetars. There are some observational clues for the existence of a magnetism-powered 
particle wind. The total particle luminosity is estimated to be ~ lO^^ergs"^, comparable to 
their persistent X-ray luminosities. Such a particle wind will amplify the star's rotational 
energy loss rate. The consequent dipole magnetic field is about ten times smaller than that 
of magnetic dipole braking, if the particle fiow is strongly collimated at the star surface. 
In the wind braking scenario, magnetars are neutron stars with strong multipole field. For 
some sources, a strong dipole field may be no longer necessary. Wind braking of magnetars 
may help us to explain some challenging observations of magnetars. 

A magnetism-powered pulsar wind nebula and a braking index smaller than three are the 
two predictions of the wind braking models. Future studies will tell us whether magnetars 
are wind braking or magnetic dipole braking. 
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